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We confirm that in order to account for the recent D0 resuh of large like-sign dimuon charge 
asymmetry, a considerable large new physics effect in is required in addition to a large CP 
violating phase in Ba — Bs mixing. In the Randall-Sundrum model of warped geometry, where 
the fermion fields reside in the bulk, new sources of fiavor and CP violation are obtained. We 
analyze the like-sign dimuon asymmetry in this class of model as an example of the desired new 
physics. We show that the wrong-charge asymmetry, , which is related to the dimuon asymmetry, 
is significantly altered compared to the standard model value. However, experimental limits from 
AMs, APs as well as K mixing and electroweak corrections constrain it to be greater than a a away 
from its experimental average value. This model cannot fully account for the D0 anomaly due to 
its inability to generate a sufficient new contribution to the width difference T\2y even though the 
model can generate large contribution to the mass difference Mf2. 

PACS numbers: 



I. INTRODUCTION 

The B factories, BaBar and Belle, have firmly established the CKM mechanism as the leading order 
contributor to CP violating phenomena in the quark sector. New physics (NP) effects can add to the leading 
order term producing deviations from the standard model (SM) predictions. These deviations are expected 
to be more pronounced in rare FCNC processes, as they are suppressed in the SM. The Belle experiment 
is scheduled for an upgrade [ll which will result in very precise results in B decays. LHCb is ready to take 
data and is expected to make many important measurements in h quark decays. These measurements may 
reveal the presence of new physics. 

In recent years, there have been several measurements of B decays which differ from the predictions of 
the SM by ~ 2a. For example, in i? — -kK, the SM has difficulty accounting for all the experimental 
measurements . The measured indirect (mixing-induced) CP asymmetry in several 6 — >■ s penguin decays 
is not found to be identical to that in B'^ J/^Ks counter to expectations of the SM and could be 
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providing hints for new physics [3]. The large transverse polarization in some penguin dominated decays 
to light vector particles, like B (f>K* are also somewhat difficult to understand in the SM where 
naively one expects the transverse polarization amplitudes to be suppressed. A further effect has recently 
been seen in the lepton sector: in the exclusive decay B'^ — > K* fi^ the forward-backward asymmetry 
has been found to deviate somewhat from the predictions of the SM 0, 0] ■ Although this disagreement is 
not statistically significant, the Belle experiment itself claims this measurement shows a clear hint of physics 
beyond the SM [8|. There are also other measurements like the branching ratio oi B ^ tv measured at 
Belle which appear to be in conflict with SM expectation 

Most discrepancies reported above have appeared in 5 — > s transitions and so it is obvious that measure- 
ments in Bg mixing will be crucial in testing the SM and finding evidence of new physics. In the SM, B^ — B^ 
mixing is generated at loop level and is suppressed. Many new physics models can contribute to B^ - B° 
mixing and can cause measurable deviations from the SM. There are already measurements in the B'^ — 
system where the mass difference AMg and the width difference AFg between the two mass eigenstates have 
been measured. Two other measurements in the Bs system have generated enormous interest as they do 
not appear to agree with the SM predictions. The first measurement is the phase of -B^ — B^ mixing which 
can be measured via indirect CP violation in Bg ~> J/4'4'- The CDF 10] and D0 11] Collaborations have 
measured indirect CP violation in Bg — >■ J/ifjcj). The experiments measured S^^ — —2(3s, and found j^l 

A=0.4lt°:i« or 1.161^^3. (1) 

This disagrees with the SM prediction 

/3f^ = 0.019 ±0.001. (2) 

Implications of this measurement for NP models have been analyzed jl^, [l^ . 

The second measurement was made recently in the i?" — B^ system when the D0 Collaboration measured 
the like-sign dimuon charge asymmetry with 6.1 fb^^ of data T4]. The following result was reported: 

A'li ^ -(9.57±2.51(stat)±1.46(syst)) X 10"^ (3) 

The like-sign dimuon charge asymmetry for semileptonic decays of b hadrons produced in pp collision is 
defined as 

where and iV^ are the number of events containing two b hadrons that decay semileptonically into 

/Lt+jU+X and ix~ ii~ X, respectively. The semileptonic decays of both Bd and Bg can contribute to A^;. The 



relation between A^i and the "wrong-charge" asymmetries af; and a^; is given by [15 1 



All = ( , ^Ti + ^1% Vsi. (5) 

where Zq = 1/(1 — y^) — 1/(1 + x^) (q = d, s). Here fd and fg denote the production fraction of Bd and Bg, 
and the quantities Xq and yq are given as, 

Xq = yq EE (6) 

where AMq and AF^ are the mass and width differences in the 5° — B^ system. The semileptonic wrong- 
charge asymmetry a^; is defined as 

, ^ TjBq -> /i+X) - T{Bq ^ ^i-X) 

r{Bq^ fi+X)+T{Bq^ f,-Xy ^' 



3 



Using the known valuesL = 0.323 ±0.037, fs ^ 0.118 ±0.015, Xd = 0.774 ± 0.008, yd « 0, = 26.2 ±0.5 



and ys = 0.0046 ± 0.027 



14 1, , one can rewrite Eq. ([5|) as 

(0.506 ± 0.043)af, + (0.49 ± 0.043)a^, 



A" — 
^si — 



The SM predictions for the charge asymmetries are [1J| 

af, = (-4.811:0) X 10-^ al, = (2.1 ± 0.6) x 10"^ 
The SM result for can be obtained using Eqs. ([8|) and (|9]) as 



(-2.3 



+0 X 10- 



(8) 



(9) 



(10) 



which is about 3.2 a away from the value in Eq. ([H]). 

The SM prediction of the charge symmetry a^; in Eq. ([9]) is consistent with the observed value a^^ = 
0.0047 ± 0.0046 [1,0, 

within errors. In order to obtain D0 measurement for yl^^ in Eq. ^ using the 
measured af; , the value of the charge symmetry a^; needs to be [l4 1 



all = -(14.6 ±7.5) X 10"^ (11) 

This value is much larger than its SM prediction in Eq. ([9]). The D0 direct measurement of a^^ = — (1.7 ± 
9.1{stat)tli{sys))l0-^ [ll], is consistent with the SM value in Eq. 
extracted by combining the D0 and CDF [3] measurements as [l^ 



An average value for a^j can be 



slJavg 



-(12.7±5.0) X 10- 



(12) 



This average value of is about 2.5(t away from its SM value in Eq. [HI A confirmation of this deviation 
would be an unambiguous evidence for new physics, and alread y interpretations of this result in terms of 
NP have been performed in various extensions of the SM 
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In this work we consider the warped extra dimension Randall- Sundrum (RS) model 23|. This model was 
proposed to solve the hierarchy problem in the SM and in this framework some of the flavor puzzles in the 
SM can be addres sed in the split fermion scenario with the fermions located at different points in the extra 
dimension 



24 



26| 



In this paper we work out the contribution to the parameters ^iid -^^12 ^^"3 -^s ~ system for the 
general case of NP with operators that are of the vector and/or axial vector types. The general formula that 
we derive can be used for several extensions of the SM. Taking the RS model as an example for new physics 
we use our general expressions to compute the contribution to Ffj and M^2- 

The paper is organized as follows. In the first section, we present an overview of the phenomenology of 
the -0° — system including constraints on NP with present measurements. In the second section, we 
present the general expression for Ffj and M(2 for general new physics containing vector and /or axial 
vector operators. FCNC effects in the RS model with split fermions are discussed in the next section. The 
subsequent sections contain our numerical results and conclusions. 



II. MODEL INDEPENDENT ANALYSIS OF B° - B2 MIXING 

In this section we will briefly review the phenomenology of the — B^ system for q = s,d. The formalism 
for B mixing is well known but we will review it here for completeness and study the constraints on NP 
imposed by measurements in this system. 
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The and states can mix in the presence of weak interactions. The resuhing mass eigenstates can 
differ in their masses and hfetimes. In the Bq — B^ system, the time evolution of the general state is governed 
by the Schrodinger equation 



dt \Bq{t)^ 



— T~Ln 



Ait) ' 



(13) 



where the Hamiltonian is given in terms of the 2x2 Hermitian mass {Mg) and the decay width (Fg 



matrices 



^'^11 2-*^ 11 ^"12 2^12 
/\/r<2* i-nq* njq i p9 
^^^^12 - 2^ 12 ^^^11 ~ 2^ 11 



(14) 



The mass eigenstates are the eigenvectors of Hq. The eigenvectors with the lightest and heaviest mass 
eigenvalues can be written as 



\B!^) = P \Bq) + q \B°q) , \B^) = p \Bq) - q , 



(15) 



with + |(7p = 1. The masses and widths of these mass eigenstates are 



Mf^ = Aff,±i?e[^(Mf2-^r?,)(MS-^r?;)], 



rf^ = rt,T2irn[MM!,-'-ri,)[M!;-'-rii)]. 



(16) 



One can now construct the following observables 



Mr, 



Af f + M^ 
—2 2- = M?- 



r - _2_ 

11' 9 ^ 



^ 11' 



A/f - M^ = 2Re\ 



AfS-^r?;)], 



:Aff2"^ny]. 



(17) 



The mass difference, the width difference and the parameters in the eigenvectors expression in Eq. (|15p can 
be written as 



iAMqr~\{ATqf 
AMqATq 



p/q 



4(|A/' 



9 |2 
I2I 



^iRe[Mf^Tlli 
'AAf, + fAr/ 



(18) 



One usually defines the two dimensionless quantities 



AAf„ 



xb„ 



ATq 

2r„ 



(19) 



Measurements indicate ub^ ^ 0(10 ^) while xb^ ^ 1- These results model independently imply 



AT, « AAf,. 



(20) 
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Thus to a good approximation 



AM, = Ar,.-i^^^ = 2|r„|co.,„ (i)^ = -^(l-i.„.|i^,), (21) 



where 



Mi2 \Mi2[^,^^ , / Ml 

'12 



0,=ar5(-|^). (22) 



The semileptonic wrong-charge asymmetry a^j is defined as 



|ri2l _ Ar, ^^^^^^ ^2^^ 



IM12I AM, 

The off-diagonal element Af of the matrix M is related to the dispersive part of the Ai? = 2 transition 
amplitude 

Ml. - ^(^.l^e^/n^")' (24) 
where tob, is the Bq meson mass. In the presence of NP contributions to it can be written as 

= Aff2^*'-l-M«2'^^ = Af«2^*'i?X/e"^'s (25) 

where 

Rli = 11 + rl,^'^' I - v'l + 2^MCOs5^, + (rl,)2, = arg\l + ^.e^^M], (26) 

with r^^e'^w = Mf^^ /M^^^'^ ■ The mass difference in Eq. (|21l) is modified as 

AAf, = AA/|f (27) 

The experimental result of AA/^ is AAf^ = 0.507 ± 0.004 ps~^, which is consistent with the SM expectation 
with i?^ ~ 1. This imposes a stringent constraint on any NP contribution to the B^ — B^ mixing. One may 
estimate the SM contribution to AAfs through the ratio AAff *^/AAf^*^, in order to minimize the hadronic 
uncertainties. For quark mixing angle 7 = 67 °, one finds AATf*^ ~ 19 ps"\ which is in agreement with 
the latest measurements by CDF [27| and D0 [3]: 

AA4 = 17.77 ±0.10(stat.)±0.07(syst.), 

AAfs = 18.53 ±0.93(stat.)±0.30(syst.). (28) 

Fig. [T] shows the allowed ranges for — r^,^ where we have neglected the SM phase in Af ^2 • The green 
scatter points satisfy the combined result of CDF and D0 in Eq. (l28l) within the Icr limit. The phase 5%^ 
is varied in the range [0, 27r], and it is not constrained below r\,j < 0.4. Also, one can conclude that R%j is 
limited by: 

0-7<Rm<IA. (29) 

The off-diagonal element Ffj of the matrix F is related to the absorptive part of the transition amplitude 
from Bq to Bq. can be written as 

r?2 = ^(^^l"^!^?)' (30) 
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FIG. 1: The allowed ranges for Sf.j — r%j from the combined result of CDF and D0 in Eq. (|28p are shown. 



where the transition operator T is defined by 



In the presence of NP contributions to F^j? it can be written as 



^ 12 



-^q,SM 



where 



\l + ry^r\= i + 24cos6l + {rl )2, 



cl>l^arg[l + ry% 



with r^e'^r = T^2^/Tfj'. The width difference in Eq.[2T]is modified as 



Ar, = 2|r?f*^|i?^cos(0|^, + 01, 



(31) 
(32) 
(33) 

(34) 



The decay width difference AFg has been measured independently. The angular analysis of Bs — >■ J/ip4> 
gives 0, 0, 

AF, = ±(0.154l0:0fo) ps-i , (35) 

to be compared with the SM prediction [iH 

AFf'^ = (0.096 ± 0.039) ps"^ . (36) 

The measurement of AF^, in principle, constrains NP contributions to Ff2- Note that the present theory 
predictions are consistent with experimental measurements though it should be kept in mind that theory 
predictions for AF^ can contain hadronic uncertainties and constraints on NP from this measurement is not 
that strong. 

In the presence of NP contributions to both and F^2, the charge asymmetry in Eq. ([23|) can be 
rewritten using Eq. ((25)) and Eq. (|32]) as 



'SM + <?^M - <?^r) 



(37) 



If one neglects the NP effects to the AB — 1 effective Hamiltonian, i.e., Ffj = Fjj^*^ (i?f. ~ 0,4>f, ~ 0), 
then using Eq. ([57)) , the charged asymmetry a^i is given by 



1 



12 



), 



(38) 
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where we have neglected the SM phase 



'SM 



= 213, 



SAI 



-0.0047 ± 0.0046, one finds that in order to account for the D0 resuhs, a 

all = (-14.6 ± 7.5) X 10~^ 



As indicated earher, using the experimental value for 

must be given by 



This implies that 



sm ( 



-(2.9± 1.5)|i?; 



(39) 



(40) 



Using the fact that |i?Af, | — 1, one finds that sin^M^ ^ 1, which is unphysical. Therefore, one concludes 
that new physics contribution to is necessar y to explain the observed CP charge asymmetry a^;. This 
conclusion was also discussed earlier in Refs. lB"22|. The measured AAf^, AF^, and {ali)(avg) in Eqs. (pS)) . 
([55]). and ([T^ . respectively can be used to determine model independently the ranges for the NP quantities 



involved in Eq. p5|) and ([32]) . In Fig. [2] we show the possible ranges for S' 



M 



' M 



sm 



\)~R 



and sin ( 



~-^) — Rm,- The green scatter points satisfy the measured AMg, AFg, and {ali)(^avg) within 
the 1(7 limit. The allowed ranges for Sfj — r\j and R\i remain the same as shown earlier. The phase 5^ is 
varied in the range [0, 27r], and it is highly constrained below < 0.4. Also, one can see that the allowed 
ranges for R^ > 1.4. These results indicate that a considerably large NP effect in Ff2 is required to address 
the observed charge asymmetry a^;. We note that similar plots and similar conclusions can be found in 
recent literature [22 1. 




1.0 1.5 2.0 2.5 3.0 3.5 



FIG. 2: The allowed ranges for 5l[ - r%f, (5f - rf, 
points satisfy the measured AMs, AFs, and (isOcaug) 



sm ( 



_Rp , and sin ( 



within the Icr limit. 



Rk 



The scatter 



III. EFFECTS OF NEW PHYSICS: GENERAL CASE 
A. NP contribution to the decay width 

In this section we present the SM and NP calculations of the off-diagonal element Ffj of the matrix Fj 
defined in Eq. (|30l). The results for Ffg including general NP are new to the best of our knowledge. We will 
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only consider vector/axial vector operators in the NP Hamiltonian. In the B'^ to transition amphtude 





FIG. 3; A graphical representation of Feynman diagrams contribute to the leading order is shown. 

in Eq. pO|) . the dominant contribution comes from intermediate cc states (see Fig. [3|)at the tree level. 
Neglecting CKM-suppressed terms one can write the effective Hamiltonian for AS = 1 transitions in the 
SM as 



32| 



AB=1 
eff 



^v*v 



6 

E 

ri=l 



h.c) , 



(41) 



where the operators are 



Ql = (biCj)v-A{CjSi)v-A, 

Qs ^ (biSi)v-A{qjqj)v-A, 
Q5 = ibiSi)v-A{q]qj)v+A, 



Q2 = {biCi)v~AiCjSj)v^A, 

Qi = {hiSj)v-A{q]qi)v-A, 
Qe — {biSj)v-A{qjqi)v+A- 



(42) 



Here the i and j denote color indices and a summation over these color indices is implied. The notation 
ifif2)v±A means fij^{l ± 75)/2 while q — u,d,s,c,b. Qi, Q2 are color suppressed and color allowed 
tree-level operators, respectively, and Qs-.-Qe are QCD penguin operators. The Wilson coefficients of these 
operators are denote d by CV..C^respectively. is determined up to next-to- leading order in both A/mf, 



and as{mb) in [3l|, [32|, [33|, [3j|, and [35|. In the rest of this section, we summarize the leading order 



results of Ffj in the l/nib expansion and introduce the effects of new physics to Ffj. We point out that in 
our numerical analysis, we used the update SM value of AFf*^ in Eq. (|36)) . 

In the SM, Ff2 for the tree-level operators Qi and Q2 is obtained by computing the related matrix elements 
in Eq. (|3T|) . The corresponding Feynman diagram is illustrated in Fig. |3l Ff2 can be written, to leading 
order in the l/m?, expansion, as j33 |. [s^ 



ri2^'' = -T^^^0^)(^c,Vcsr{Go{z) <Q-> +Gso <Qs->), (43) 

where z = ml /ml and the functions Gq{z) and Gso are given by 

Go(z) = VT^[{1-z)Ga + ^{1-^z)Gb), 
Gsoiz) = VT^^{1 + 2z){Ga-Cb). (44) 



The combinations of Wilson coefficients Ga and Gb are given by 

Ga = (tVcC? + 2C1C2), Gb^GI, 



(45) 
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where Nc ~ 3. The two AS ~ 2 operators Q_ and Qs- hi Eq. (|43l) are given by 

Q_ (6,7^^(1 - 75)s»)(&,7m(1 - 75)sj), Qs- = (&,(1 - - 75)5,). (46) 

The matrix elements < Q- > and < Qs- > can be found using a vacuum insertion method as (ssf 

< Qs_ > = Qs_ = -l^ljlj^^B^, (47) 
where _Bi and B2 are bag parameters, and the quantity is defined by 



ml. 



(48) 



^ (m{,(mb) + ms(m6))2' 
The contributions of the all six operators Qi{i — 1..6) in Eq. to Ffj are given by 

ri2^*' = - ^J(|^|^) (k;^c.)^(go(^) < q- > +g;jo(^) < qs- > ), (49) 

where the functions G'q{z) and G'gQ{z) are given by [s 



2 

G's^{z) - VT^(l + 2z)(C;-C^). (50) 
The combinations of Wilson coefficients and G'^ are given by 

C^i - (iVc(Ci+C3)2 + 2(Ci+C3)(C2 + C4) + 2C5C6 + iVcC2), 

= ((C2+C4)2 + C|), 
G'c = 2(iVc(Ci+C3)C5 + (Ci + C73)C6 + (C2 + C4)C5 + (C2 + C4)C6). (51) 

Setting Ci=3..6 = in Eq. one recovers Eq. The effects of Ci=3..6 on V2i^^ are small. 

Next we consider the contributions to Ffj from NP operators involving b — )■ ccs transitions. Note that 
these operators can also contribute to the lifetimes of the h hadrons and we will consider constraints on 
these operators from lifetime measurements in our numerical analysis. In general NP can also contribute to 
h — )■ qqs transitions with q ^ u,d, s. However, we note that in the SM these transitions are suppressed and 
measurements in decays like B — >■ Kn, B — > (pK etc. already constrain these NP contributions. In other 
words, an NP contribution to 5 — )■ qqs with q = u,d, s cannot significantly affect and hence we ignore 
these contributions. The b — > ccs transitions, on the other hand, are tree level in the SM and measurements 
of decays with the underlying b — > ccs transitions, such as Bd — > DgD, J/tpKs etc. might still allow for NP 
in 6 — !■ ccs transitions that are not small. In this section, we only present model independent results, and 
for the numerics in the following section we consider the RS model with split ferniions. The RS model with 
split fermions can also generate NP contributions to 6 — >■ st~^t~ which are not that well constrained from 
experiments. However, these transitions are generated by the exchange of KK electroweak bosons and hence 
they are much smaller than the b — )■ ccs transitions that are generated by KK gluon exchange. 

The AB = 1 effective weak Hamiltonian for NP operators is written as 

= {^^llQll + >^'llQ'll + ^lrQlr + >^'lrQ'lr 
+^rrQrr + ^'rrQ'rr + AfliQ RL + ^RlQ RL^ I 

(52) 
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where 

QlL = (biSi)v-A{CjCj)v-A, Q'lL = ibiSj)v-AiCjCi)v~A, 

QlR = {biS.i)v-A{CjCj)v+A, Q'lR = {biSj)v -A{cjCi)v +A, 

Qrr ~ {biSi)v+A{cjCj)v+A, Q'rr. = ibiSj)v+AicjCi)v+A, 

Qrl = (biSi)v+AicjCj)v~A, Q'rl = {biSj)v+A{cjCi)v-A- (53) 

The eight coupUngs Xab, and X'j^g (A, B = L, R) are in general complex and can be determined from 
specific models. Thus, the total effective Hamiltonian can be written as 

nj _ njSM , njNP 

Heff = ^i^cbVcs) ({Cl + C3 + C3")Ql + (C2 + C4 + Cf ^)Q2 + {C5 + Ct")Q5 + (^6 + Cf'')Qe 



+C^''Qrr + CrQkfl + Q'^^QflL + C«^Q'flL j , (54) 
where the new Wilson coefficients (A, B = L, R) are 

^iAA _ ^^AA ^AA _ '^^AA 



^AB _ V^^AB ^AB ^ ^^'aB 



Ci"" - 7^^, Ct^^^^^ {A^B). (55) 



The NP contributions to T21 are obtained by computing the related matrix elements in Eq. (|3T|) to the 
leading order in the l/rrib expansion using the MS scheme. The NP contributions contain both pure NP 
and SM-NP interference terms. In general the latter dominate NP contributions due to the large SM Wilson 
coefficients. We obtain NP contributions to T21 as 

^12 = ^12 +^12 +^12 : (56) 

where T^i^^ , T^'^^ , and Y^'™^^ contain the contributions from LL and LR, RR and RL, and all possible 
type of operators, respectively. They can be expressed in terms of the matrix elements of eight AB = 2 
operators 

Q- = (biSi)v~A(bjSj)v-A, Qs- = (biSi)s^p(bjSj)s~p, 
Q+ = {hsi)v+A(bjSj)v+A, Qs+ = (biSt)s+p(bjSj)s+p, 
Qt = (biSi)v-A(bjSj)v+A, Qs^ = (biSi)s-pibjSj)s+p, 

Q± = {btSi)v+A{bjSj)v-A, Qs± = {biSi)s+p{bj^j)s-P, (57) 
where S± P = 1 ±75. The explicit forms of Vl^^^ , T^^'^^, and p^^™*^ are 

^12" = - ^2^g^^) (K;K.)^(G-(z) < g_ > +G'^o(.) < Qs- > 



1f - -T^:;^0^)(^ctVcsr{Go{z) < Q+ > +Gso{z) < Qs+ : 

if" = -y^^^(K;K.)^(Gi(z) + G2(z) + G3(z) + G4(z)), (58) 
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where the G functions are 
G^'(z) 



Gsoiz) 
Gi(z) 



VT^((1 - z)C'; + i(l - 4z)C'^ + 3zC'i), 



VT^( (1 - z)Ca + - 4z)Cb + 3zCc 



^3 



+ cf«)) < > 



VT^(1 + 2z)(Ca-Cb), 

2VT^(Gi + C3 + Cl'^) [((1 - zHNcC^'^ + Q^^) + 3z{NcCi 
+(1 + 2z)(7VcC5^^ + Ce^^) < QsT > " . 
G2(z) = 2^/T^(G2 + Ci + Gf^) [(szGg^^ + (1 - z)G5^^) < > +(1 + 2^)Gf^ < Qs^ > 

-i(l + 2z)G6«^ < Q± > -2(3zGf « + (1 - z)G6«^) < Q5± > ] , 
Gsiz) = 2VT^{C5 + Gl-^) [((1 - z)(7VcG3«« + Gf ^) + SziNcC^"^ + Gg^^)) < > 

+(1 + 2z) (iVcG3«« + Gf ^) < Qs^ > , 
G4(z) = 2VT^{Ce + Gg^^) [((1 - z)G3^^ + SzGg^^) < > 

+ (1 + 2z)G3«« < Qs^ > +3z(-Gf « < Q± > +2Ge«^ < Qs± >) • 



(59) 



The combinations of Wilson coefficients are 



C'c = 

Ca = 

Cb = 

Cc = 



Nc[{C^n' + {Ci^'T) + 2G3^^G4^^ + 2G5^«G6^'^^ + 2G3^^(7Vc(Gi + G3) + G2 + G4) 
+2G4^^(Gi + G3) + 2C^"iNcC^ + Ge) + 2G6^«G. 
(G4^^)2 + {Cl^^f + 2G4^^(G2 + G4) + 2G6^^G6) 

2G3^^(7VcG|'^ + Gl-^) + 2G4^^(G|'^ + Gg^^) + 2G3^^(iVcG|'^ + C^^) 
+2G4^^(G5 + Ge) + 2G|'^(iVc(Gi + G3) + G2 + G4) + 2G|'^(Gi + G2 + G3 + G4) 
NcUci'^f + {Ci'^f) + 2G3^^G4^^' 



RR\2 



2{NcC^^C^'^ + C^'^C^'^ + Gi^^Gs^^ + G4^^G6«^) 
The matrix elements < Q >— {Bs\Q\Bs) of the operators in Eq. (|57|) are given by [3 

<Q-> = <Q+>= ^mljl^Bi, 

5 

<Qs-> = < Qs+ >= -i;m%jly^B2, 



(60) 



< QsT > = < Qs± 2m%Ji^y^Bi 

4 



< 



> = < Q± >= -^mljlj^^B^, 



(61) 



where B^'s are the bag parameters and their numerical values for 5° — i?^ system in the Af S'-NDR scheme 
at ruf) — 4.6 Gev can be found in 37| . 



12 



B. NP contribution to the mass parameter MI2 

In this section we present the SM and NP calculations of the off-diagonal element Mfj of the matrix Mq 
defined in Eq. (12^ . The effective Hamiltonian for A_B = 2 transition that generates B^-B^^ mixing in the 
SM can be written as 

^AB=2.SM ^ 3_M'^{Vt,V:,fllB^S,{xt){bs)v^A{hs)v^A, (62) 

where t]b^ — 0.551 [s^ is the QCD correction, and the loop function £'0(2:4) {xt — m^/M^) is given by [4 



Axt-llx^t+x^t ixllnjxt) 



The SM contribution to MI2 can be obtained using Eq. (l24l) as 

^^12^*' = -^M2,(y,fcF;j2^B.5o(xt)mBj| Bi. (64) 



Next we consider NP effects in the AB — 2 transition. The effective Hamiltonian for this transition is 
written as 



(65) 



where (A, B = L, R) are NP couplings. The operators Q^, Q+, Q^, and Q± are given in Eq. ([57)1 and 
their color suppressed counterpart operators can be written as 

Q' = {biSj)v~A(bjSi)v-A, Q'+ = (bi3j)v+A(b]Si)v+A, 

Q'^ = (biS.j)v-A(bySi)v+A, Q'± = (biS.j)v+A(bySi)v-A- (66) 

Using Eq. (|24p and applying the Fierz transformation, one can obtain the contribution of these NP operators 
to as 



12, NP 



^ {{^'ll + 5ll) <Q-> +{S'rr + 6rr) < g_ > +{Slr + Srl) < g=F > 



-2{S'^j, + 6'j,j^)<Qs^>). (67) 
For certain classes of models including the RS model with split fermions the following relations hold: 

Sll = -1/3S'^L, Srr = -1/3^^^, Slr = -I/SS'lr, Srl = -1/3(5^^, S'^r = S'^^- (68) 
One can then obtain 



4 



2/„, „, \„ /I 



3 

where the bag parameters are defined in Eq. (jSl 



IV. NEW SOURCE OF FLAVOR AND CP VIOLATION IN THE RS MODEL 



The warped extra dimension model has been proposed as a solution of the hierarchy problem. In the 
original RS model, the SM fields are localized to one of the boundaries and gravity is allowed to propagate 
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in the bulk. However, it was realized that scenarios with SM gauge bosons and fermions in the bulk may 
lead to a new geometrical interpretation for the hierarchy of quark and lepton masses. The Higgs field has 
to be confined to the TeV brane in order to obtain the observable masses of the W and Z gauge bosons. 
We will consider the scenario of Ref . 23] , based on the metric 



(70) 



where cr(y) — K\y\ and k ^ Mp is the curvature scale determined by the negative cosmological constant in 
the five dimensional bulk. The fermion fields reside in the bulk of this nonfactorizable geometry and can be 
decomposed as 



(71) 



Here Tc is the radius of the compactified fifth dimension on an orbifold S1/Z2 so that the bulk is a slice of 
^^5*5 space between two four dimensional boundaries. The left-handed zero mode wave function is given by 



24| 



/f(c/„,l/) 



-c/o<T(y) 



where c/^ = mf^/ k, and m^^ is the bulk mass term. Using the orthonormal condition: 

1 



dye'^(^)/f (c/„,y)/f (c/„,y) = l, 



one finds that iVo is given by 



/ 7rKrc(l — 2c/^) 
The right-handed zero mode wave function can be obtained from 



(72) 



(73) 



(74) 



(75) 



The tower of fermion KK excited states is not relevant to our discussion here. Note that fermions with 
c/ > 1/2 are localized near the Planck brane at y = and fermions with c/ < 1/2 are localized near the 
TeV brane sX y = nrc- 



The massless gauge fields that propagate in this curved background can be decomposed as [24 1 

CO 

A,,(x,y) = -=^4")(a;)/(")(y), 



(76) 



with fl^\y) is dimensionless. The n-th (n > 0) mode function is given as 



gj(a) 



An) 



Jl 



!i4_g<T(y) 



+ 5^(m^'Vi 



An) 



'A _^a(v) 



where Ji and Yi are the J- and Y-type Bessel functions of order one and 

(„) _ Jo(to^"Vk) 



bA{m) 



(77) 



(78) 
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The coupling of the gauge KK modes to the fermion at the vertex <}2(r)1l{r)9''^'' given by 

.(5) 



5' 



(27rrc)3/2 



e'^^^^ icu , , {y)dy. 



(79) 



The nonuniversal parameters c/^ lead to nonuniversal couplings to the KK state of the gluon. In the basis 
of mass eigenstates we have the following flavor dependent couplings: 



fjTU.d{n)\ /-w7-u,d'^\ (n) / \ /-^ruA \ 

\'^L{B) l-^P - i^L(i?,)i"7 5L(fl)(.C/J (,l/L(fl)^7/3: 



where the g^^jj) is given by 
(1) / \ ( 1 - 2^j„ 



(80) 



(81) 



The tree-level relation between the 5D and 4D QCD couplings 175 and g is g = g^jsflwFc. However, at the 
one loop level the relation between these two couplings also depends on the value of the brane-localized 
kinetic terms for the bulk gauge fields, and there can be significant corrections to the tree level relation. A 



detailed discussion on this topic can be found in Ref. 41 1. Finally, the unitary matrices V^'^'^j diagonalize 



the up/down quark mass matrix M^^'^ = (w/v^)^^^'^, which is given in this model as 
The dimensionless parameters lap are defined as 

lap = >^al^^ 



(82) 



(83) 



(5) 

where A^,^ are the Yukawa couplings which are free parameters to be fixed by the observable masses and 
mixing. In Table HI we present an example of the c/ parameters that leads to the correct quark masses and 
mixing with A^^-* ~ C'(l). The corresponding first-KK gluon coupling constants are given by 



<7i'^(cQj 


= -.199, 




-.198, 


5l^(cq3) 


= 1.496, 




= -.191, 


5l^(cd2) = 


-.191, 


5l^(cd3) 


= -0.198, 


5k (ct/J 


= -.199, 


9l {cu2 ) = 


-.195, 


5l (cc/3) 


= 3.38. 



(84) 

We note that the value of in Table. |T]is larger than for the 1st /2nd generation. Such a choice of c's is not 
consistent with anarchic 5D Yukawa couplings discussed in Ref [26]. We also point out that for simplicity, 
we have used here the tree-level relation between the 5D and 4D QCD couplings. As can be seen from these 
values, (/^^■' is of order one. In fact, this general conclusion can be obtained with any other values of c/. One 
can easily show that for c/ > 1/2 the coupling 5'-'^' approaches zero, while for c/ < 1/2 one finds 5(1) < 4. 
We expect the general features about t;'^^ to remain true even when the one loop matching of the 5D and 
4D QCD couplings is used. 

In the RS model, the effective Hamiltonian for AB = 2 transition can be generated at tree-level via the 
exchange of a KK gluon, as shown in Fig. ID The AB = 2 effective Hamiltonian is given by 



42 



n 



AB=2,KK 
eff 



E 

A,B 



4m 



32 



KK 



(85) 



15 



9 



in) 



FIG. 4: The gluon KK contribution to B° - B° mixing is shown. 



with A, B = L, R, and Pl{r) ~ 1/2(1 T 75)- The contribution to MI2 in the RS model can be obtained by 
comparing Eq. (|85p with Eq. ()65|) . and the coupUngs can be expressed as, 



jd(i) 



16m 



KK 



^LR - te::^ — ^'^L myUR j32, - ^^^2^^ )32[u^ ) 



(86) 



Thus, Eq. (IMl) reduces to 

■KK 1 



12m 



2 "^bJI^ 

KK 



((t/f + (C/h^'^)^2)Si + iui^%2K^%2{lr^^B, - 3r^i34)] • (87) 

This result agrees with those in [s^. If M'l^^ is dominated by the first term, then the ratio tm^ in Eq- 
(PS)) can be obtained for the RS model as 

/2415x2/([/f^' 
V mjfK ^ ^ Vti 



^KK 

Ms 



,KK 



M- 



s,SM 
12 



2^2 



rd(l)^ 



Therefore, r^f - 0(1) requires ttikk 2.4 TeV if we assume {Uf^')z2 ^ VthVts- If (t^'''^^)32 =i C(l), 
the experimental limits on AMs^ implies that vtikk ^ lOTeV, which imposes stringent constraint on the 
associated compactification scale. 

The AS = 1 effective Hamiltonian for the h — ^ sec transition in the RS Model can be written as 

1 '-r'^(l)^„^^r^"(l)^„J^^;.,,MD.„.^^^.,, i( 



e// 



4^2 -(C^r'032(t/B^'^)22 (6a^PASj)(c,7Mi^BQ) - - (6.7^^^^^, ) (c, 7m^bc, ; 



''KK 



(89) 



The contribution of the RS model to can be obtained using Eq. ([56l) with the couplings 



LL 



S'r 



LR 



S'rr = 



^'rl = 



1 


(L/,^W 


16m|.^ 


1 


(L/,^W 


16m|.^ 


1 




16m|.^ 


1 




16m|.^ 



(90) 



The corresponding Wilson coefficients in the RS Model can be obtained from Eq. 
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NUMERICAL RESULTS 



The numerical inputs for the parameters in the SM 16[ and RS model are summarized in Table. ID The 
values of the bag parameters in the M^-NDR scheme can be found in 37[, and the decay constant of Bg 
is from [s^. The relevant CKM matrix elements are obtained from the CKMfit collaboration [43j. The SM 
Wilson coefficients for the quark level b — )• sq'q transition at next-to-leading order in the NDR scheme are 
obtained from 



44| 



In the RS model, the matrix elements M}^'^^ and Yl^'^^ depend on the four couplings {U'^^-^)32, and 



(f/2(fl))22 [see Eq. m] 



Writing these couplings explicitly yield 





)32 


^ ^L(32) *^L(22) 


[9^^ 




-9^L\CQ^)]- 








)32 






(CD2) 


(1)/ M 


-^<33)<32)[5i'^(cA,) 


(1),' M 




22 


~ *^L(21) *^i(12) 


\9'P 


(cQi) 


~9^l\^Q.)]~ 


^<22)^L(22)[ffi'^(cQj 


-9l\cq^)]^ 




22 


~ ^R(2l) ^B.(12 


\9'^ 


(cc/i) 


-9r\cu3)] H 


^ yR(22) ^fl(22) [5^?^^ (cc/2 ) 


-9r\cuJ\, 



(91) 

where the unitarity of the V^^^ is used. The bulk parameters c/^ specify the position of the fermion's 
localized wavefunction in the bulk. The specific choice of c's are model dependent and they should generate 
the 4-d Yukawa hierarchy of the quarks as well as the CKM mixing in the left-handed sector. In our 
numerical analysis, we do not consider any specific values for c's. We scan the allowed values for the 
couplings (J7^|^^-))32, and (C/^|]j^-))22 after imposing the constraints from the experimental measurements of 



AMs and AF^. Assuming the matrices to be CKM- like and using g'^^^ < 4 one notes from Eq. (|9T|) . 
that the coupling (C/^|^-,)32 is constrained to (C/^|^^-,)32 ^ 0.1, while (?72(ij))22 is of order one. Therefore, we 

vary the four KK couplings as : |(C^^(^^))32| ^ 0.1, \{U2''^^)22\ ^ 1, and |(J7]j'"^')22| ^ 3, and their phases in 
the range [0, 2tt]. Note that this assumption oi Ul and bein g CK M-like is consistent with the scenario of 
anarchic 5D Yukawa couplings that has been considered in Ref. 26| . Our choice for the right-handed mixing 
angles are not quite consistent with the scenario of anarchic 5D Yukawa couplings but our assumption is 
that with proper tuning of the 5D Yukawa couplings we can have the mixing angles in the range considered 
in this analysis. 

In the fit, AMs and AFj are constrained by their experimental results within Icr errors. All SM input 
parameters are uniformly varied within their errors and the SM Wilson coefficients are evaluated at fif, = 4.2 
GeV. The bag parameters are kept at their central values and niKK is varied in the range [1.2, 10.0] TeV. 
Note that strong constraints on to ^-y are obtained from measurements in K mixing [45[ but there are still 
regions of parameter space ( 38|, |4l|, |46| ) which allow values for uikk considered in this work. Note that the 
lower end of the KK scale which is scanned here (a couple of TeV) would require some amount of tuning in 
order to satisfy the constraint from ex (even with Higgs in the bulk)[33|- Similarly, the KK scale of 1.2 to 



2 TeV used in the scan is generally disfavored by electroweak precision tests 2l|, |47| . Such a low KK scale 



mig ht be allowed if the 1st and 2nd generation fermions are all chosen to have a (very) close-to-fiat profile 
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but then one loses the explanation of fermion mass hierarchies based solely on profiles. 



We observe that AMs(exp) constrains \{uf}^^)32\ to < 10 ^, and allows r 



Note that similar conclusion about r^f^^ was also reached in Ref. 



,KK 
M 



\M. 
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_ . The fit results allow the ratio 
^s,KK _ ips^/f/f^ps^S-A/i ^ \Q% for values of ttikk < 1-5 TeV and this ratio falls quite quickly as mxK is 
increased beyond 1.5 TeV. As indicated earlier, NP in & — )• ccs transitions will also contribute to the lifetimes 
of h hadrons. We expect the corrections to the total widths of the b hadrons to be also < 10% which cannot 
be detected experimentally given the hadronic uncertainties in calculating the total widths 48| . The lifetime 
ratio of the Bg meson to the Bd meson, has a tiny theoretical uncertainty in the SM [48|. New physics 
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TABLE I: Numerical values of the theoretical quantities used in the numerical analysis are shown. 



0.00 



-0.01 



Numerical values for the input parameters. 


rriBg = o.oDt) (jeV 


JBs = 26o(y.O) MeV 


mf,(mt) — 4.19(+U.l8)(-0.UD) GeV 


i3l=0.87 


mc(mc) = 1.27(+0.07)(-0.09) GeV 


-62=0.84 


ms(i GeV) = 0.01 GeV 


-63=0.91 


ms{mi,) =U.0o4 GeV 


-64=1.16 


tb^ ~ 1.425 (0.041) ps 


-65 = 1.75 


It/ I nr\yiioo 

|\/ci, = 0.04128 


CQi = 0.72 


IKcsl = 0.97J42 


= U.b 


\Vtb\ = 0.999141 


CQg = 0.35 


l^tsl = 0.04054 


— 0.63 


Ci(mt,)=-0.1903 


C!72 = 0-30 


C2{mb)= 1.081 


CU3 = 0.10 


C3(m6)= 0.0137 


= 0.57 


C4(mt,)=-0.036 


= 0.57 


C5(m6)= 0.009 


CD3 = 0.60 


C6(m6)=-0.042 






^•-^1.0 -0.5 0.0 

sin[0^-4] 



0.006 
0.004 
0.002 
0.000 
-0.002 
-0.004 
-0.006 




4 6 8 
»iKK[TeV] 
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FIG. 5: The a|;- sin ((ji^j^ — ) (left panel) and a^j- rriKK (right panel) correlation plots in the RS model. The gray 
band indicates the la experimental allowed ranges for a^i (avg). The green line indicates the SM prediction for a^;. 



will contribute equally to the Bg and Ba total widths, in the leading order in the heavy quark expansion, 
thereby largely canceling in the lifetime ratio. Our naive estimate is, NP in b — >■ ccs transition can correct 
this lifetime ratio as — « ^^(1 + x)^ with x = { 1 — pSA/ where F^" j are the SM Bg d widths. Using 

rp^g- < 10 %, = IzbO.Ol IS, we get x < 0.1% which is consistent with experimental measurements 161. 
In Fig. [5] is shown the a^; - sin — (pf,) and a^i-viKK correlation plots in the RS model. The gray band 
indicates the la experimental allowed ranges for a*; (avg) while the green line indicates the SM prediction 
for all- As one can see from these figures, the fit results allow a^i ~ —0.00498, which is 1.54(t away from its 
experimental average value in Ea. ([T2)) . The value of the corresponding sin (0m — </>f) —0.76. Hence this 
model cannot fully account for the experimental results and this is due to the fact that this model cannot 
generate enough correction to the width difference AFfa. Note that, as rriKK is increased the suppression 
to a^i due to rnKj^ is partially compensated by larger mixing angles which are within the considered ranges 
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in this work and are consistent with experimental measurements of AMs and AFj. However, after a certain 
point the suppression due to ttikk cannot be compensated and a^; decreases with increasing uikk mass. 



VI. CONCLUSIONS 



In the past few years several measurements in rare B decays involving & — s transitions have been 
somewhat difficult to understand in the SM. This has put the focus on measurements in the B'^ — system. 
The measurements of the phase in B^ — B'^ mixing and more recently the like-sign dimuon asymmetry have 
generated an enormous interest in the community as these measurements indicate possible deviations from the 
SM predictions. In this work, we presented calculations for general new physics corrections to the parameters 
in — B'^ mixing. Taking the Randall-Sundrum model as an example of new physics we calculated the 
wrong-charge asymmetry, a^j, as well as other parameters in B'^ — B^ mixing. Our calculations indicate 
that while the RS model can cause deviations from the SM predictions for the wrong-charge asymmetry, 
it cannot explain the present experimental average value within the la range. This is due its inability to 
generate sufficient new contribution to the width difference AFf 2 , even though the model can generate large 
contribution to the mass difference AMfn. 
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